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CMB lensing extraction and primordial non-Gaussianity 
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The next generation of CMB experiments should get a better handle on cosmological parameters 
by mapping the weak lensing deflection field, which is separable from primary anisotropies thanks 
to the non-Gaussianity induced by lensing. However, the generation of perturbations in the Early 
Universe also produces a level of non-Gaussianity which is known to be small, but can contribute 
to the anisotropy trispectrum at the same level as lensing. In this work, we study whether the 
primordial non-Gaussianity can mask the lensing statistics. We concentrate only on the "tempera- 
ture quadratic estimator" of lensing, which will be nearly optimal for the Planck satellite, and work 
in the fiat-sky approximation. We find that primordial non-Gaussianity contaminates the deflec- 
tion field estimator by roughly (0.1/jvl)% at large angular scales, which represents at most a 10% 
contribution, not sufficient to threaten lensing extraction, but enough to be taken into account. 

PACS numbers: 98.80.Cq 



Introduction - Cosmic Microwave Background (CMB) 
anisotropies are of considerable interest for cosmology be- 
cause after cleaning the observed temperature and polar- 
ization maps from various foregrounds, one obtains a pic- 
ture of cosmological perturbations on our last-scattering 
surface. The power spectra of primary anisotropies are 
related to various cosmological parameters, and depend 
on the physical evolution mainly before the time of de- 
coupling (and also, more weakly, on its recent evolution, 
through the integrated Sachs- Wolfe effect and the angu- 
lar diameter-redshift relation). 

It was realized recently that CMB anisotropies encode 
even more cosmological information than expected, be- 
cause it should be possible in a near future to measure 
the deflection field caused by the weak lensing of CMB 
photons by the large scale structure of the neighboring 
universe at typical redshifts ^ ~ 3 P, 0|- The power 
spectrum of the deflection field encodes some informa- 
tion concerning structure formation mainly in the linear 
or quasi-linear regime, and is therefore extremely useful 
for measuring parameters like the total neutrino mass 
or the dark energy equation-of-state, which mildly affect 
the primary anisotropy So, the next generation of 
CMB experiments could output for free a Large Scale 
Structure (LSS) power spectrum, without suffering like 
galaxy redshift surveys from the systematics induced by 
mass-to-light bias and by strong non-linear corrections 
on small scales at z < 0.2. 

There are several methods on the market for extract- 
ing the deflection map 0, |E IE I3Lj all based on the non- 
Gaussianity induced by lensing Q . These methods start 
from the assumption that both the primary anisotropies 
and the deflection field are Gaussian; they also assume 
that the noise present in the temperature and polariza- 
tion maps is Gaussian and uncorrelated with the signal. 

None of these assumptions is exactly true. Amblard 
et al. II] already estimated to which extent the lensing 
extraction will be biased, first, by the non-Gaussianity 



of the lensing potential caused by the non-linear growth 
of matter perturbations on small scales, and second, by 
the imperfect cleaning of the CMB maps from the ki- 
netic Sunyaev-Zel'dovich effect, which also has a black- 
body spectrum, induces non-Gaussianity, and features 
spatial correlations with many of the structures respon- 
sible for the lensing. Both effects were found to be rele- 
vant (i.e., to induce a significant bias in the estimators). 
However, they are small enough to preserve the validity 
of the method. 

The purpose of this work is to relax the first of the pre- 
vious assumptions, and to consider realistic situations, in 
which one cannot avoid a small level of non-Gaussianity 
to be produced in the Early Universe. Non-Gaussianity 
emerges as a key observable to discriminate among com- 
peting scenarios for the generation of cosmological per- 
turbations and is one of the primary targets of present 
and future CMB satellite missions Indeed, despite the 
simplicity of the infiationary paradigm [lol |. the mecha- 
nism by which cosmological curvature perturbations are 
generated is not yet established. In the standard slow- 
roll scenario associated to one-single field models of infia- 
tion, the observed density perturbations are due to fluc- 
tuations of the inflaton field itself when it slowly rolls 
down along its potential. In the curvaton mechanism |T]| 
the final curvature perturbation is produced from an ini- 
tial isocurvature perturbation associated with the quan- 
tum fiuctuations of a light scalar field (other than the 
infiaton), the curvaton, whose energy density is negligi- 
ble during infiation. Recently, other mechanisms for the 
generation of cosmological perturbations have been pro- 
posed, the inhomogeneous reheating scenario the 
ghost inflationary scenario and the D-cceleration 

scenario just to mention a few. Single-field slow-roll 
infiationary model inflation itself produces a negligible 
amount of non-Gaussianity, and the dominant contribu- 
tion comes from the evolution of the ubiquitous second- 
order perturbations after inflation. However, alternative 
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models for the generation of perturbations might pro- 
duce much stronger primordial non-Gaussianity. There- 
fore, non-Gaussianity in the CMB maps from primordial 
fluctuations could mask the non-Gaussianity from lens- 
ing distortions. We will see later that if one expands the 
power spectrum of the lensing estimator in powers of the 
gravitational potential $ ~ 10~^, the contribution from 
primordial non-Gaussianity appears at the same order as 
the lensing power spectrum itself. Therefore, a precise 
computation is needed in order to understand whether 
the primordial non-Gaussianity could affect lensing ex- 
traction. 

Lensing extraction with quadratic estimators - Weak 
lensing induces a deflection field d, i.e., a mapping be- 
tween the direction of a given point on the last scatter- 
ing surface and the direction in which we observe it. At 
leading order this deflection field can be written as 
the gradient a lensing potential, d = V(/). 

In the limit of Gaussian primordial fluctuations, the 
unlensed anisotropics obey Gaussian statistics, and in 
the flat-sky approximation their two-dimensional Fourier 
modes are fully described by the power spectra Cf^ where 
a and b belong to the {T, E, E\ basis. Weak lensing cor- 
relates the lensed multipoles 0, ^3 according to 



ensures that d"** is an unbiased estimator of the lensing 
potential: 



(d'^^(L))cMB = ^L0(L) = d(L) 



(5) 



Note that, so far, the coefficients g'^^(\i,\2) are still ar- 
bitrary. From the observed temperature and polariza- 
tion maps, one could compute each mode of d"** and ob- 
tain various estimates of the deflection modes, precise up 
to cosmic variance and experimental errors. In order to 
quantify the total error, it is necessary to compute the 
power spectra of the quadratic estimators 



d°''*(L)d°''(L)) = (27r)25(L - L')C; 



dd{ab) 
L 



(6) 



where the average is now taken over both CMB and LSS 
realizations, since 0(L) is also a stochastic quantity. In 
this definition, the power spectra are written with a su- 
perscript dd(ab) in order to be distinguished from the 
actual power spectrum of the true deflection field. These 
spectra feature the four-point correlation function of the 
observed (lensed) Fourier modes (a(li)6(l2)a(l3)6(l4)), 
which should be expanded at order two in (/)(L) in or- 
der to catch the leading non-Gaussian contribution. 

The four-point correlation functions are composed as 
usual of a connected and an unconnected piece. The con- 
nected piece is by definition a function of the power spec- 
(a(l)6(l'))cMB = (27r)2^(l + l')C','''' + /"''(l, + (1) tra C?*' in which we now include all sources of variance: 



where the average holds over different realizations (or 
different Hubble patches) of a given cosmological model 
with fixed primordial spectrum and background evolu- 
tion (i.e. fixed cosmological parameters). In this aver- 
age, the lensing potential is also kept fixed by conven- 
tion, which makes sense because the CMB anisotropics 
and LSS that we observe in our past light-cone are sta- 
tistically independent, at least as long as we neglect the 
integrated Sachs- Wolfe effect. The above function 
is defined in and takes a simple form in the case 
ab = TT: 



cosmic variance, lensing contribution and experimental 
noise. The unconnected piece is a function of the same 
spectra plus the deflection spectrum Cf"*, and as usual it 
can be decomposed in three terms corresponding to the 
different pairings of the four indices ^^: {lijh), ihjh) 
or (lijl^), (hjh) or (^1,^4), (l^jh)- The first term leads 
to considerable simplifications when it is plugged into 
the expression of the quadratic estimator power spec- 
trum, and the result is simply Cf"^, as one would expect 
naively from squaring Eq. ||SJ). The other terms lead to 
more complicated expressions that we will write as two 
noise terms: 



/^^(i,i') = cr(i+i')-i+cr(i+i')-i'. (2) 

Our study will be based on the quadratic estimator 
method of Hu & Okamoto (which is equivalent in 

terms of precision to the alternative iterative estimator 
method of Hirata & Seljak as long as CMB exper- 
iments will make noise-dominated measurements of the 
B-mode, i.e., at least for the next decade). By invert- 
ing Eq. (^), one builds a quadratic combination of the 
temperature and polarization observed Fourier modes 



^dd{ab) _ ^dd 



ab 
cJL 



(7) 



which represent respectively the contribution from the 
connected piece and from the two non-trivial terms of the 
unconnected piece [Hill (later, we will give the exact 
expressions in the case aa = TT). In order to get an 
efficient estimator, we should adopt the set of coefficients 
(7°''(li,l2) which minimize the noise terms. It is actually 
much easier to minimize the connected term only, which 
leads to the simple results 



d"''(L) 



a(li)&(l2)5"''(li,l2) (3) 



where I2 = L — li, and in which the normalization con- 
dition 



/tt(1i,12)5"''(1i,12) 



(4) 



<?""(ll,l2) 



/aa(l,10 



and [N, 



aa 
cJL 



J^aa 



(8) 



for a ~ b (for a ^ b see @). With such a choice, the 
unconnected piece contribution Nf' can be shown to be 
smaller than Al'^, but not completely negligible [T^ . 

The various estimators d"'' can be constructed for 
each pair of modes, except for the pair BB, because 
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the spectrum is dominated by lensing at least on 

small scales, which invalidates the present method. So, 
the quadratic estimator technique would not be optimal 
for long-term CMB experiments with cosmic-variance- 
dominated measurement of the B mode 0, ■ For an 
experiment of given sensitivity, the five other estimators 
can be combined into a final minimum variance estima- 
tor, which gives the best possible estimate of the deflec- 
tion field by weighing each estimator accordingly to its 
noise level. The sensitivity of the Planck satellite ^3 is 
slightly above the threshold for successful lensing extrac- 
tion, but only at intermediate angular scales, and with 
essentially all the signal coming from the d^-'" estimator. 
The following generation of experiments - such as the 
CMBpol or Inflation probe project [23| - should obtain 
the lowest noise level from the d^-^ estimator 0. 

We summarized here the quadratic estimator method, 
which assumes that both the primary anisotropics and 
the lensing potential are Gaussian. We will now study 
the impact of primordial non-Gaussianity. This program 
is numerically cumbersome, and we will only concentrate 
on the d"^-^ estimator, which is the only relevant one for 
Planck, while sticking to the flat-sky approximation in 
which numerical computations arc much quicker. 

Contributions from primordial non-Gaussianity - The 
two-dimensional Fourier modes of both temperature 
anisotropics and the lensing potential can be related to 
the stochastic three-dimensional modes of the primor- 
dial gravitational potential <&(k), multiplied by a transfer 
function which accounts for its time-evolution. The fact 
of writing a unique stochastic function can bring some 
confusion, because the modes $(k) which appear in the 
CMB and lensing expressions represent fluctuations at 
very different redshifts: the first ones before decoupling, 
the second at z ~ 3, i.e. in the neighboring universe. So, 
as long as we neglect the integrated Sachs- Wolfe effect, 
it is convenient to introduce two statistically indepen- 
dent functions $'-^^^(k) and $'"^^(k), sharing the same 
statistical properties, but sourcing respectively a(l) and 
</)(!). 

The true non-Gaussian potential ^E"]^^ (x ~ CMB or 
LSS) can be expanded in real space in powers of a Gaus- 
sian potential (f>^. In Fourier space and at order three. 



$^Jk) = <i>^(k) + $5(k) + $^(k) 



(9) 



d3 



with 

cl>f (k) *(p,)$f *(p2)<i>f (pi+P2+k) 



P $f(k+p)ci>f*(p)-(27r)35(k)<i>i 



and 



(27r)3(27r)3 



potential. They are proportional to the dimensionless 
parameters /nl and (/nl, respectively. The theoreti- 
cally predicted parameter /nl appears as a kernel in 
Fourier space, rather than a constant, in most of the 
scenarios for the generation of the cosmological pertur- 
bations, while theoretical predictions for the parameter 
(7nl arc still lacking This gives rise to an angular 

modulation of the quadratic non-linearity, which might 
be used to search for specific signatures of inflationary 
non-Gaussianity in the CMB. In this paper, however, we 
restrict ourselves to the simplest case and assume /nl 
and gNL as mere phenomenological multiplicative con- 
stants. Under this assumption, the WMAP team has 
measured the bispectrum to obtain the tightest limit to 
date, -58 < /nl < 134 (95%) 113 . On the other hand, 
no observational bound has been set on ^nl from the ob- 
served trispectrum. However, one can simply notice from 
Eq. Q that the small parameter (/nl*!*) contributes at 
the same order as (^(?nl so, by comparing with the 
/nl bound, it is likely that values of order y^g^L ~ 100 
are still allowed by the data. 

As far as lower bounds are concerned, one should keep 
in mind that although single-field slow-roll inflation it- 
self produces a negligible amount of non-Gaussianity, 
the dominant contribution comes from the evolution of 
the ubiquitous second-order perturbations after inflation. 
This effect must exist regardless of the inflationary model, 
setting the minimum level of non-Gaussianity in the cos- 
mological perturbations at order /nl^5jvl^1- 

In order to evaluate the impact of these extra contri- 
butions on the power spectrum (^^'^t"''')^ should first 
recompute the four-point functions (a™'6"^'^a™'''6™''), 
working as before at order six in the gravitational poten- 
tial. Non-zero contributions can arise only from terms 
with an even number of <I>2'^^(k) and (E>^^^(k) factors. 
The standard calculation of the previous section included 
terms in which either two multipoles were lensed at order 
one in 0(1), or one mifltipole was lensed at order two (the 
later terms contributes only to the connected piece). In 
addition, we should now consider terms in which: 

A. the four multipoles are unlensed, but two of them 
include the term $5^^, 

B. the four multipoles are unlensed, one of them in- 
cludes the term ^g'^^, 

C. one of the four multipoles is lensed at order one in 



(/)(1), which includes the term 



alss 



We have parametrized the primordial non-Gaussianity 
by a quadratic and a cubic term in the gravitational 



The last term C vanishes because (p^f^ has zero av- 
erage. So, at leading order, lensing and primordial non- 
Gaussianity effects are completely separable, and we sim- 
ply need to add corrections from the primordial non- 
Gaussianity trispectrum, which is given in Okamoto & 
Hu (who computed it in the Sachs- Wolfe approxima- 
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FIG. 1: 



(1), for the case of Planck (left) and 

CMBpol (right). The thick (red) curve shows the variance of the signal Cf"^ . Other curves represent (from top to bottom at 



Various contributions to the variance of a single mode of the estimator d 



I ^ 2000) the noise variance from the connected part of the four-point correlation function, from the lensing trispectrum, and 
from the terms B, A2 and Al in the primordial non-Gaussianity trispectrum (displayed for /jv^ = g]^^ = 100) . 



tion). First, one needs to compute the integrals 



Fi{ri,r2) 
ai{r) 



- / k^dkP^{k)ji{kn)ji{kr2) , (10) 



k^dkAi{k)ji{kr) , 
k^dkP^{k)Ai{k)ji{kr) 



(11) 
(12) 



where A; (fc) is the radiation transfer function for the tem- 
perature, normalized to $ = 1 in the early universe. We 
checked that our functions ai{r) and /3i{r) (computed 
with a slightmy modified version of CMBFAST |2^) per- 
fectly agree with those of . Each of the A and B-type 
trispectra are composed of three parts, which can be sim- 
ply expressed in terms of the intermediate quantities 



[(^h{ri)Ph{ri) + c.p.] [ai^{r2)f3u{r2) + c.p.; 



\hl2 



2/2 r^dr [ai,ir)l3i,{r)f3i,ir)PUr) + c.p. 



where c.p. means circular permutation of the indices 
li. The final expression for the power spectrum of the 
quadratic estimator for any mode L of modulus L, in- 
cluding the three trispectra induced by lensing and pri- 
mordial non-Gaussianity, reads 



C 



dd{TT) 



Ar 





fdHi d% 




1 (27r)2(27r)2 



g^^ih,h)9^^il[X2) 



2|i,-i;r2C|lf_,,|/^^(u,~i;)/^^(i2,-i^) 



3(pb);?i? 

L-, Lr, 



(13) 



with I2 = L— li and I'j = L— I'j^. The variance of the signal 
Cf''' comes from the first term of the lensing trispectrum, 
while the noise variance A^^ comes from the connected 
part of the four-point correlation function. We can iden- 
tify all the other sources of noise inside the integral. The 
third line of Eq. (|13|) contains the two other terms of the 
lensing trispectrum, which give equal contributions (as 
can be shown after a change of variable). The last line 
contains the three terms of the A-type trispectrum (the 
last two contribute equally) and the three terms of the 
-B-type trispectrum (three equal contributions). Let us 
call Al, A2 and B the terms of the last line. 

The numerical integration of Eq. (|13() could be ex- 
ceedingly long, due to the four-dimensional integral in 
Fourier space (which converges only when the upper 
bound of integration is taken around the Silk-damping 
cut-off, Zinax ^ 2500) and of the two-dimensional inte- 
gral in r space. Fortunately, the second integral can be 
limited to values of r close to the comoving distance to 
the last-scattering surface, far from which r'^ai{r) is neg- 
ligible. Here, we will integrate r = tq — t from 0.03t» 
to 1.3t, with a conservative step dr — O.Olr, (where r 
stands for conformal time, tq is evaluated today, and r, 
is the conformal time at recombination, defined as the 
peak of the visibility function) . 

The integration in Z-space can be sped up by taking 
into account various symmetries, and also by noting that 
for the terms A\ and B the integral over li and l'^ are 
separable. Then, the integration effectively reduces to a 
two-dimensional one. For these terms, we integrate first 
on li, then on r\ and r2, and obtain the full result in few 
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minutes. The term A2 and the lensing term require a 
real four-dimensional integral and more CPU time. For- 
tunately, the quantities to integrate are very smooth in 
Z-space, and one can considerably reduce the computing 
time by choosing a large step without loosing precision. 
We checked that Al ~ 30 is by far sufficient. 

Results and conclusions ~ We take a fiducial ACDM 
model with f7fc = 0.05, fic = 0.25, f^A^O.TO, /i = 0.65, no 
reionization and a scale-invariant primordial spectrum 
P,^{k) = 6.204 X 10~^^fc~'^. For instrumental noise, we 
consider the cases of Planck HFI (three channels) and 
of the CMBpol project, with a sensitivity described by 
the same parameters as in We show in Fig. ^ the 

various contributions to the estimator power spectrum, 
as computed from Eq. ifT^ . Note that we are plotting 
the variance of a single mode 1, and not the error on the 
reconstructed deflection power spectrum, which can be 
lowered by combining all modes of given wavenumber I 
and by binning the data (this is why Planck is likely to 
make a reasonable detection of the deflection power spec- 
trum at intermediate Ts 0|, although the noise variance 
is slightly larger than the signal variance in Fig.^. The 
contributions from the primordial non-Gaussianity terms 
A and B scale respectively like and gjvL, and here 
they are shown for /^^ = \J9nl = 100. 

The contamination from primordial non-Gaussianity 
appears to arise mainly at low from the Al-type term. 
On these scales it would be necessary to perform an exact 



all-sky computation in order to make a precise prediction. 
However, the error caused by the flat-sky approximation 
even at low I is usually small ^3 • 

The noise induced by primordial non-Gaussianity is 
responsible for roughly (0.1/7vl)% of the amplitude of 
the estimator d^-'"(l) in the range 2 < ^ < 10: so, around 
10% for the largest possible value of /atl, and around 
0.1% for standard slow-roll inflationary models. In the 
range 100 < I < 1000, the contribution is roughly of 
order {IQ-^ fNL)% from the A2 term, and (0.01^5]^)% 
from the B term. 

If in the near future /nl appears to be large, it will be 
measured independently using the three-point correlation 
function. It should then be possible to subtract to some 
extent the bias induced by primordial non-Gaussianity 
when reconstructing the power spectrum C'lf, as sug- 
gested in ll,!!^ for other sources of bias. 

We conclude that primordial non-Gaussianity should 
be taken into account if it is as large as to saturate the 
present upper bounds, but that in no case it will represent 
a dangerous issue for lensing extraction. 
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